Dynamical phase transition of a ID transport process including death 
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Motivated by biological aspects related to fungus growth, we consider the competition of growth 
and corrosion. We study a modification of the totally asymmetric exclusion process, including the 
probabilities of injection a and death of the last particle 5. The system presents a phase transition at 
Sc{a), where the average position of the last particle (L) grows as Vi. For 5 > 5c, a non-equilibrium 
stationary state exists while for 5 < Sc the asymptotic state presents a low density and max current 
phases. We discuss the scaling of the density and current profiles for parallel and sequential updates. 
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Fungi are eukaryotic organisms that include micro- 
organisms such as yeasts and molds, as well as the fa- 
miliar mushrooms. The growth of a fungus is formed by 
the combination of the apical growth and the branching 
process which leads to the development of a mycelium. 
Most of them grow as hyphae, which are a cylindrical 
thread-like structures of 5 — 10 um in diameter and up 
to several centimeters in length [l| . The apical growth is 
a quasi one-dimensional process that extends the hypha 
by transport of material from the seed to the front tip. 
At this later position, enzymes are released into the en- 
vironment, where the new wall material is synthesized. 
The rate of extension, in a favorable environment, can be 
extremely rapid, up to 40 micrometers per minute. Many 
bio-physical models describing the growth of fungal 
colonies and/or single hypha have been studied. In this 
context, efforts focused on the non-equilibrium proper- 
ties of a modification of the totally asymmetric exclusion 
process (TASEP), by consideri ng a distinct dynamics of 
one of the two boundary sites [THllj. 
An important aspect not taken into account yet, is that 
fungi have the ability to grow in a wide range of habitats, 
including extreme environments 12|-|l4j and survive in- 
tense UV/cosmic radiation during space travel fl5!| . Since 
the wall of the tip is usually structurally weak 1] , in such 
a situation the extension rate can be slowed down and as 
the hypha is progressively aging, it may break down or be 
broken by other organisms (l| . Our analysis is focused on 
the theoretical description of the above-mentioned grow- 
ing process in competition with a corrosive environment. 
The theoretical model used is based on a simple modifi- 
cation of the TASEP and captures the general behavior 
induced by the two competing processes. 
Proposed in 1968 to study the motion of ribosomes along 
mRNA [H, nil, numerous modifications of the TASEP 
were introduced including multilane systems and multi 
species transport [18, 22], Langmuir dynamics 23-25], 
extended par ticles 0, [irj as well as systems with finite 
resources 2^, 2^. In a general framework, such models 



are considered as toy models of transport phenomena in 



FIG. 1: Sketch of the microscopic rules of the system. A 
particle is injected on the first site (if empty) with probability 
a while the last particle dies with probability 5. In the bulk, 
particles are jumping forward with probability one if the site 
is empty. 



order to better understands physics far from equilibrium. 

In the first part of this letter we define the model 
and the system parameters. After a careful definition 
of the parallel update dynamics we obtain, in the 
stationary state, the exact expression of the generating 
function and discuss the scaling of the density profile. 
In the last part, the results for the sequential update 
are presented. We are focusing our attention on the 
density and current profiles used to characterize the 
phase diagram. The system exhibits a dynamical 
phase transition between a region for which the size 
of the chain reaches a stationary value {{L)st) and 
a region in which the length is increasing with a 
constant velocity. At the transition line, the system 
exhibits a diffusive behavior {L){t) cx y/t. Finally we 
summarize our results and discuss areas for future work. 

Our model is a dynamical extended exclusion process 
(DEEP) described on a semi-infinite one-dimensional lat- 
tice. The dynamics is defined by the probability a of 
injection of a particle (on the first lattice site, if empty) 
and the probability 6 of death of the last particle. In the 
bulk, the particles are hopping to the right with prob- 
ability p — 1, while the last particle (if not destroyed) 
jumps from site L to L + 1 with probability 1 — ^ (see fig- 
ure [1]). As a result of the microscopic rules, the length of 
the active chain (L) (defined by the rightmost occupied 
site) presents a dynamical phase transition. Extensions 
of this model can be considered including partially asym- 
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metric diffusion and/or absorption/desorption processes 
of particles, so called Langrauir dynamics. 

Considering first the case of the parallel update dy- 
namics, we define the site occupation of the site i at time 
t by r]{i,t). The position of the last particle, at time 
t, is denoted by L such that all sites i > L are empty 
(r]{i,t) = 0). In one timestep, the last particle is up- 
dated first, it can jump forward or die. Then all the 
particles in the bulk are simultaneously moved forward 
and finally a particle is injected on the first lattice site 
with probability a. Therefore the probability to have n 
empty sites between two particles is given by (1 — a)". 
Knowing the configuration at time t, given by the set of 
ri{i,t) {i G N*), the dynamical rules are: 

rj{L + l,t + 1) ~ 1 with probability 1 — S 
77(1 < i < L,t + 1) ^ ri{i — 1, t) with probability 1 
r/(l,t + l) = l with probability a. 

After one timestep, the probability to find the last par- 
ticle at the same position L is given by the product aS. 
This situation appears if we kill the last particle which 
is immediately followed by another one. By extension, 
knowing that ri{L, t) — 1, the probability to find, at time 
t -I- 1, the last particle in L — n is given by 5{\ — a)^a. A 
transition is expected between two regions separated by 
the line 5c{o)- For 5 > 5c, the last particle is killed at 
a sufficiently high frequency, that the growing process is 
effectively stopped. In this case a stationary state exists 
and we have {L)st < 00. On the other hand, below the 
critical line the growth is only slowed down by the death 
of the last particle. The growing rate is finite, there is no 
stationary state and {L){t) oc t. In the later situation we 
will analyze the asymptotic density and current profiles. 
Starting in the region 5 > 5^, as we approach the critical 
line, the average {L)st diverges. We can then neglect the 
effect of the left boundary and write the master equation 
for the probability P{L,t) to find the last particle at site 
L at time t. This describes a biased random walk with 
probability I — 6 to jump to the right while the walker can 
jump to the left, over n sites, with probability Sa{l — a)"'. 
We can write explicitly 

P{L,t+l) = {l-S)P{L-l,t) + SaP{L,t) 

00 

+ fo^(l-a)"P(L-|-fc,i). (1) 

k=l 

The time evolution of the average size is given by (L) (t -|- 
1) — {L){t) = 1 — S/a and leads to the equation for the 
critical line Sc = a. For 6 > Sc, &n exact expression of the 
growing rate is given by the velocity v = (L) jt = X — bja. 
At the critical point, the equation for the second moment 
leads to a diffusive law ((i^) oc t) such that (L) = D^^^/t. 
The coefficient -D|| has been obtained from simulation 
data and evaluated to be (see figure [3|) 



The discrete time version of equation ([Ij , with the appro- 
priate boundary equations (for L = 0, 1) can be cast into 
a single equation for its generating function defined by 
t) = z^P{L,t). Analysis of the long time limit 
shows for S < Sc that (L) — >• (1 — S/a)t is recovered 
and leads, for 6 > Sc, to the exact stationary expression 



l-a-z{l-S) 



(3) 



Using {L)st = dzG* {z)\z=i we obtain the expression of 
the stationary length {L)st — a(l ~ a)/ [S{S — a)]. The 
order parameter {L)~^^ vanishes linearly as we approach 
the critical line {{L)~^ cx e with e = 5 ~ Sc ^ t). Finally 
the stationary probability P*{L) (for L 7^ 0) is 



1-^ 

1 — a 



(4) 



with P*(0) — 1 — a/S. The characteristic length A = 
[ln((l — a)/(l — diverges at the critical point. 

The exact expression of the density profile, in the sta- 
tionary state, is given by the sum p*{n) = aX^L ^ 
n)P{L), where 9 is the Heaviside function. This leads to 
p*{n) = {S/a){N)stP{n), with {N)st = E„>iP*H such 
that the average density p = {N)st/{L)st is equal to a. 
At the critical point, the scaling structure of the proba- 
bility P{L,t) is suggested by the diffusion law (L) cx \/t 
and leads to P{L, t) ~ t^^^"^ f{L/ \/t). For large time, the 
asymptotic form is f{u) cx exp(— u^/7rZ?jj). For 5 < Sc, 
the distribution is gaussian, centered in (L) {t) such that 
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P{L, t) cx — = exp 



L-{L)(t) 



(5) 



By integration we obtain the density profile given by 



p{x,t) 



1-erf 



x-{L){t) 



(6) 



(2) 



where erf(a;) is the error function. The expression 
of the critical line could have been obtained from 
the equation j = j^, where j = a is the particle 
current arriving at the interface while js = 5 \s the par- 
ticle current leaving the system due to the death process. 

For sequential update dynamics, at each timestep, 
M{t) + 1 updates are realized by choosing at random 
among the N{t) particles (present on the lattice at time 
t) and the first lattice site. Then on average, during 
one timestep each particle is updated once. When the 
rightmost particle is updated it can be removed from the 
lattice or moved forward. If a particle of the bulk is se- 
lected it moves systematically forward with respect to 
the exclusion process. Finally when the first lattice site 
is updated, a particle is injected with probability a if 
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the site is empty. The analysis of the system behavior 
reveals a more complex phase diagram showed figure [2] 
In the region S < Sc, the chain size (L) is diverging and, 
for large time, the system properties should be similar 
to an infinite TASEP with effective /3e// = 1- The crit- 
ical line is then obtained when the current of particles 
arriving at the interface (j) is equal to the current of 
particles leaving the system {js = 6). Since, for 6 < 6c, 
the system is expected to behave like the TASEP another 
transition line should exist for a — 1/2, between the low 
density and max current phases. The current of parti- 
cles is given by j — a{l — a) for a < 1/2 and j = 1/4 
for a > 1/2. We claim that the exact expression of the 
critical line is given by Sc = a{l — a) for a < 1/2 and 
Sc ~ 1/4 for a > 1/2. This is confirmed by the Monte 
Carlo simulations (fig. [5]). In order to characterize these 
two phases, we are focusing our attention on the density 
and current profiles. 

It is seen (figure |4]) that the quadratic relation which 
usually expresses the current as a function of the local 
density of particles, j = p{l — p), does not hold close to 
the interface {x ~ (L) ) but is recovered in the bulk. The 
correction C{x,t) defined as p{x,t){l — p{x,t)) — j{x,t) 
is shown as an insert of figure |4] for a = 1/2 and 
(5 = 1/8. For large time, C{x,t) is gaussian and rescaled 
as C{xt^^^^ — vt^^'^,t), where v is the front velocity. 
In the low density phase, the density is constant given 
hy p = a, while in the max current phase, as ob- 
served in [8|, the current and density profiles continue 
to evolve as the system size is growing. The asymptotic 
expressions of the profiles are obtained assuming x/t^ 
being the scaling variable of the density profile. Using 
the mean field relation j{i,t) = p{i,t){l — p{i + l,t)), 
one obtains in the long time and continuous space limit 
3{x,t) = p{x/V'){l - p(x/t^)) + Oit"''). Together with 
the continuity equation dtp+dxj = 0, the density is given 



by p{x/t^) 



1 — zj) which imposes z = 1. Finally, 



for large times, the bulk density and current profiles are 

p(^A) = ^ (i - f ) , jix/t)^i(^i^^y (7) 

The numerical data and analytic predictions are in per- 
fect agreement as shown in figures |3] and [S] 
In the stationary region, the probability distribution 
P*{L) as well as the density and current profiles decay 
exponentially, while on the critical line the decay of the 
density and current profiles is gaussian. Interestingly, for 
the sequential update dynamics, the coefficient of diffu- 
sion (Dj_ — {L)/^/t) presents a non monotonic evolu- 
tion as a function of a (figure [3]). In the low density 
phase, since the number of interactions between particles 
is small, the same behavior is observed for D\\ and D^. 
We numerically checked that the first correction is of or- 
der ^/a such that D± ~ D\\ + 0{^/a). However, while a 
increases, the interactions between particles occur more 
often and lead to a transition from the low density phase 




< L > (t) = D t 



Max current phase 



FIG. 2: Phase diagram of the model. The critical line be- 
tween the stationary and diverging regions is given by the 
equation 5^ = ct(l — a) for a < 1/2 and (5c = 1/4 for ct > 1/2, 
the points are obtained via Monte Carlo simulations. The 
diverging region is divided in a low density and max current 
phases. 



to the max current phase. A signature of this transition 
is found in the difference D±^ — D\\ . In the insert of figure 
131 we plot the derivative da{D± — D^^) which presents a 
maximum at the transition line a = 1/2. 

In summary, motivated by the competition between 
growth and corrosion processes, we obtained the phase 
diagram of a dynamical extended toy model inspired by 
mycology problems. For both parallel and sequential up- 
dates, the model presents a dynamical phase transition 
between a region of finite size and a diverging region. 
For the sequential update, the transition between the low 
density and max current phases leaves a signature in the 
difference D± — -D|| . It is also seen that the quadratic law 
j — p(l — p) fails close to the interface but remains valid 
in the bulk. As observed in Q, in the max current phase, 
the density and current profiles are time dependent. In 
this region, we proved that the profiles are function of 
the scaling variable x/t and independent of the bound- 
ary conditions. An identical dynamical phase transition, 
characterized by the same diffusive behavior on the tran- 
sition line, was observed in 3l| on microtubule dynamics 
for which the growth enters in competition with a detach- 
ment process. It suggests that universal features might 
emerge close to the transition line generated by the com- 
peting processes. We motivate the experimental char- 
acterization of this transition by the analysis of fungus 
growth under UV radiation used as a control parameter. 
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FIG. 3: CoefBcient of diffusion as a function of a. The red 
squares are numerical results obtained for parallel updates 
and are in agreement with the expression ((JJ - dashed line. 
The black circles, for sequential updates, present a non mono- 
tonic evolution. In the insert, the function da{D± — Dy) 
presents a signature of the transition at a = 1/2. 
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FIG. 5: Current profile as a function of the position x plotted 
for the times t = 1000, 5000, and t = 10000. The black, violet 
and green curves are obtained respectively for the parameters 
(a = 1/4, S = 3/32), {a = 1/2, S = 1/8) and (a = 3/4, 
5 — 1/8). The dashed lines are given by the expression ((7|. 




FIG. 4: Density profile as a function of the position plotted 
for the times t = 1000, 5000, and t = 10000. The black, violet 
and green curves are obtained respectively for the parameters 
[a = 1/4, S = 3/32), {a = 1/2, S = 1/8) and (q = 3/4, 
6 = 1/8). The dashed lines are given by the expression ([7|. 
In the insert, the correction C{x,t) is plotted for the times 
{t = 5000 and t = 10000) and for a = 1/2, S^l/S. 
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